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On the basis of studies of radiation-induced shape-change of reactor graphite GR- 
280, through the series of measurements of samples with different orientation of 
cutting with respect to the direction of extrusion, a conclusion is made about the 
existence of polycrystal substructural elements - domains. Domains, like graphite as 
a whole, possess the property of transverse isotropy, but have different amplitudes 
of shape-change and random orientations of the axes of axial symmetry. The model 
of graphite, constructed on the basis of the concept of domains allowed to explain 
from a unified point of view most of existing experimental data. It is shown that the 
presence of the disoriented domain structure leads to the development of radiation- 
induced stresses and to the dependence of the shape-change on the size of graphite 
samples. We derive the relation between the shape-change of finite size samples and 
the actual shape-change of macro-graphite. 
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I. INTRODUCTION 

Reactor graphites of various grades, used in nuclear engineering as moderators of neu- 
trons and constructional materials, are subjected to high dose damaging radiation. At that, 
developing in graphites radiation-induced physical, mechanical and dimensional effects sig- 
nificantly affect their operational characteristics. 

Detailed analysis of such irradiation effects as a shape-change, evolution with irradiation 
dose of elastic moduli, thermal expansion coefficient, temperature- and electrical conductiv- 
ity, radiation creep, etc., suggests an essential role of graphite morphology in all irradiation 
effects.- - - The morphology of graphite is formed by two main components - the filler and 
the binder, as well as due to the presence of an ensemble of microcracks and technological 
poresi^, and to a large extent is determined by the technology of graphite fabrication. 
The binder has a fine crystalline structure, whereas the structure of the filler is hierarchical, 
based on microcrystallites with more or less perfect hexagonal crystal lattice of two types 

- ABAB... and ABC ABC ', the fraction of which may be different.- Between the crystal- 
lite with hexagonal structure and the macro-graphite as a whole (meaning the graphite as 
a material of finite sizes) there is a succession of scaled structural levels of organization 
of crystallites, of a binder and ensembles of microcracks, which form the morphology of 
graphite.-^ 

Elastic constants of hexagonal crystallographic lattice are described by five independent 
constants en, c%2, C13, C33 and C44. 8 It is established that reactor graphites, obtained by ex- 
truding or pressing, have properties of a transversely isotropic medium 1 ' 9 , and to describe 
their elasticity they also require five independent variables, for example, E±,E\\,G\\,u_i_,v\\ 

- Young's and shear moduli and Poisson's ratio in the isotropy plane and normal to it, 
respectively. Radiation-induced shape-change of graphite has similar symmetry properties, 
allowing us to introduce the diagonal tensor of radiation- induced shape-change (see below). 

For extruded graphites the axis of axial symmetry is given by the direction of extrusion (|| 
direction), whereas directions normal to it (_L directions) define the plane of isotropy.- 1 The 
well known reactor graphite GR-28Q~ belongs to this group of graphites. In the present paper 
we discuss experimental results obtained from studies of radiation-induced shape-change of 
graphite GR-280. 
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II. DESCRIPTION OF EXPERIMENT 

Graphite samples in an amount of several hundreds obtained by longitudinal and trans- 
verse cut from the bulk material (|| and _L orientations), of diameters 6 mm, 8 mm, 12 mm, 
and of length 50 mm were exposed to irradiation at temperatures 460 ± 25°C, 550 ± 25°C, 
640 ± 25°C in the research fast-breeder reactor BOR-60 up to dose levels ~ 3 x 10 22 n/cm 2 
at energies E n > 0, 18 Mev. Measurements were carried out with the dose intervals 
~ 0.6 x 10 21 n/cm 2 . The samples were measured both in longitudinal and transverse di- 
rections, using the 4.5 mm step. When measuring the cross-section the diameters at 32 
points on the circle were recorded. 

All data are presented as functions of doses of neutrons with the spectrum of BOR-60. 
Dose conversion to the neutron spectra of thermal reactors are presented in Ref.— . The 
conversion of neutron doses BOR-60 into equivalent doses and fluxes of neutrons of thermal 
reactors requires a special consideration and will be given in a separate paper. 

As the result of mass measurements of samples, in the study of patterns of their shape- 
change the two observations were obtained, that required an additional explanation: 

• Circular cross-section of samples acquires with the irradiation dose an increasingly pro- 
nounced elliptical shape, and orientations of ellipses vary randomly for measurements 
along the sample on the scale of ~ 6 mm. 

• Results of calculation of the relative volume changes using conventional expressior*i~- 

(AV/V) C = (AL/L)" + 2 (AL/L) ± , (1) 

(where (AL/L)" and (AL/L) ± are relative length changes of the samples with parallel 
and perpendicular orientations, respectively) increasingly diverge with irradiation dose 
from results of direct measurements of the relative change of volume. 

Below we present the interpretation of the above facts in the framework of the model 
evolved by authors that takes into account hierarchical morphology of graphite. 1,1 ' Defining 
the characteristic element of a certain scale in the scaling hierarchy as radiation-induced 
afline shape-changing region, the concept of the "domain" is introduced, as a region of in- 
termediate size between the grain of graphite and the macro-graphite (of size ~ 5 mm). The 
domain has the same symmetry (but different values) of radiation-induced shape change as 
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macro-graphite, nevertheless its symmetry axis does not coincide with the axis of extrusion. 
We use this concept to develop a formalism that allows to bind together all the above facts. 

III. PROCESSING EXPERIMENTAL DATA 

In this section we describe the geometry of graphite samples before and after irradiation 
and introduce the main physical concepts to describe radiation-induced change of their 
shape. Using these concepts we find the relation between local characteristics of domain 
structure and experimentally observed radiation-induced shape-changes of graphite samples. 

Domains. Graphite is a polycrystal, and as each polycrystal it consists of small mi- 
crocrystallites. The most important difference of graphite from usual polycrystals is the 
presence of mesoscale three-dimensional objects with the size £ > 0.5 cm for brevity called 
domains, see Fig. [TJ Each domain consists of many microcrystallites, possesses the proper- 
ties of transversely isotropic medium and can be characterized by individual orientation of 
the axis of symmetry. At preparation of graphite the local orientations of anisotropy axes of 
different domains deviate from their average over the sample orientation due to technolog- 
ical reasons. The change of shape of domains under irradiation is restricted by their local 
environment, since this change is constrained by neighboring microcrystallites of different 
orientations, leading to the development of internal stresses in the bulk material, see Fig. [T] 
b. These stresses lead to the additional strong dimension effect in graphite. The knowledge 
of the dimension effect is extremely important to understand the behavior of bulk graphite 
materials, whereas most of experimental data are obtained for finite size graphite samples. 

Because of transversal isotropy of a domain its deformation can be described by two 
principal strains: along the principal axis of axial anisotropy, (Al/l),,, and perpendicular to 
it, (Al/l) ± . These strains determine the tensor of radiation-induced deformation Fq which 
is diagonal in local coordinate system, axes of which are directed along the main axes of 
local deformations: 

f l + (M/l) ± ^ 

^o= l + (A///) ± • (2) 

v l + (A///) ||y 

a) Homogeneous structure. Before irradiation the specimen has the shape of a 
cylinder of length L and diameter d. After irradiation the shape of the cylinder will be 
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FIG. 1: Graphite is a polycrystal consisting of domains (shown schematically by hexahedrons) of 
typical size £. Each domain can be characterized by the local direction of anisotropy (shown by 
small arrows), a) In the case of homogeneous structure all domains have the same orientation, b) In 
the case of heterogeneous structure these directions are random vectors with preferential orientation 
(shown by large arrow), determined by the specimen orientation. Radiation-induced shape-change 
of randomly oriented domains leads to the rise of cross-domain stresses with irradiation doze. 

changed because of radiation induced shape-change of randomly oriented domains. We 
first consider the case when orientations of axial symmetry axes of all domains of the sample 
coincide with the extrusion direction, so every sample can be characterized by the orientation 
angle to between the direction of axial symmetry and the axis of the specimen, see Figs. CD 
a and [2] a. Although the approximation of homogeneous texture can be justified only for 
small disorder in domain orientation, it is commonly used in the literature because of the 
simplicity of measurements of sample sizes. 

Relative deformations of the specimens along i = x,y,z -directions in the coordinate 
system related to the sample are described by the deformation tensor F idea i, (ALj/ / Li) ideai = 



Fideai ) — 1. In the case of homogeneous texture the deformation tensor has the form 





(3) 



where Ri (u) is the tensor of rotation by the orientation angle cu: 
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FIG. 2: Before irradiation the sample has cylindrical shape with circular cross-sections, Figure a). 
Different coordinate systems b): (x, y, z) is attached to the sample with axis z along the preferential 
direction, (x.',y',z') is related to a domain with axis z' along the direction of axial symmetry of 
the domain. After irradiation the sample takes a "crankshaft" shape with elliptical cross-sections, 
characterized by maximal d max and minimal (d m i n ) k diameters measured along main axes rotated 
by random angles (fo) k , Figure c). 

b) The structure with random domain orientations. The orientation of local 
axial axis in a cylindrical specimen can be characterized by two angles: the azimuth angle 
<Po in the cross-sectional plane (x, y) and the polar angle 9 between the local axial direction 
and the cylinder axis, see Fig. [5] b. Both these angles randomly vary in space, and their 
average values, ipo and u = 9 depend on orientation of the sample - how it was cut from 
a bulk graphite array. In the following we will consider both parallel (with average axial 
direction along the cylinder axis, u = 0°, superscript ||) and perpendicular (with average 
axial direction perpendicular to the cylinder axis, u = 90°, superscript _L) orientations. 

Macroscopic deformation of the sample is described by deformation gradient tensor F 
defined in coordinate system, related to the specimen. For the axial axis rotated by the 
angles (<po, 9) with respect to the axis of the cylinder the gradient tensor takes the form 

F = R(tp ,6)F R- 1 (tp ,6) (5) 

Here R (ip , 9) =_R 2 (fo) Ri (9) and the tensor of rotation in the plane of cross-section of the 
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sample R2 (y?o) is defined by expression 



' cos (fo — sin (p ^ 



R2 (ipo) 



simp cos^o • (6) 

V V 

The orientation of the deformation tensor F of a domain randomly varies in space on the 
length £ about the domain size. As the result of such variations of F the specimen is 
deformed after irradiation into a "crankshaft" shape, see Fig. [2] c. The bending of the 
cylindrical specimen at locations of domain junctions is described by effects of higher order 
in angle deviations and will not be studied below, although this effect is steadily fixed 
experimentally. Consider the shape of the crankshaft cylinder in more details: 

Transverse dimensions. Before irradiation the cross-section of the specimen has 
circular shape with angular dependence of the diameter do (<f) = (d cos <p) x + (dsiny>)y. 
After irradiation the shape d (99) of the cross-section is determined by the projection of 
the vector Fd {tp) on the cross-sectional plane (x, y). Calculating the square of the cross- 
sectional diameter in the deformed sample we find its dependence on azimuth angle ip: 

d 2 O) = d cos2 (<P ~ fo) + d 2 min sin 2 (ip - (p Q ) (7) 

Thus, initially circular cross-section of diameter d after irradiation takes elliptical shape, see 
Fig. [2] c. Maximal d max and minimal dmin diameters of the ellipse are given by equations: 

rfmax/rf = 1 + (Al/l) ± , (8) 

d min /d = 1 + (Al/l) ± cos 2 9 + (A///) u sin 2 9 (9) 

According to Eq. (jSJ) the maximal diameter d max is identical for all segments of the specimen, 
while minimum diameter <i m ; n randomly varies between different segments. 

Relative elongation. We denote by Lo/c the initial length of k-th domain. After 
irradiation the length of the fc-th segment of the crankshaft cylinder can be found as 
the projections of the vector F (L fcz) on the direction z of unit vector along the axis of the 
cylinder: 

L k = L ok \l + (Al/l) ± sin 2 9 + (A///) y cos 2 6\ (10) 
The total length of the specimen is found as the sum of lengths of all its segments 
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From Eqs. ( ITT]) we get the relative change of the length 



(AL/L) = (Al/l) ± sin 2 6 + (A//0|| cos 2 6 = (A///),, + e (12) 

Here (AL/L) is measured as relative elongation of the whole length of the specimen, while 
(Al/l) ± and (Al/l)^ are corresponding components of the principal strain. We use the 
notation 7TT for averaging over the sample: 

JV 

.Y 



T , - T 



k=l 

and Ek is the flattening factor of k-ih cross-section: 



sin 2 #, (13) 



^max (iin)t Ad max (Ad m [ n ) ^ , » 

^ = d = d (M) 

Volume change. The volume of the specimen can be found as the sum of volumes of 
all its segments, see Fig. [2]c: 

EN 77 7T 

^ _^ ^dmax (d m i n )^ Z/fc — ~d max d mm Z; (^^) 

Here f d max (dmin)^ is the area of fc-th cross-section. Using this equation we get the relative 
change of the volume of such "free" sample 

(AV/V) F = (Ad max /d) + (Ad~/d) + AL/L 

= 2 (Ad max /d) + AL/L - e (16) 

According to Eqs. (jSJ) and (IT21 the volume change can be expressed through the principal 
strains of the domain: 

(AV/V) F = Tr (f - l) = 2 (Al/l) ± + (A///),, (17) 



where Tr (• • ■ ) is the sum of diagonal elements. Eq. (fTTj) demonstrates the self-consistency of 
our definition of principal strains, since it reproduces the experimentally measurable volume 
change. 

Processing experimental data. In the following we will calculate the two princi- 
pal strains (AZ/Z)|| and (Al/l) ± through the change of experimentally measurable maximal 
diameter and length of the specimen: 

(AZ/Oii = AL/L - e, (Al/l) ± = Ad max /d (18) 
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where we used Eqs. (|T2|) and ([8]). In our experiments the set of diameters (di) k of the 
specimen is measured along N equally separated cross-sections k for M equidistant angles 
(fj = lixjjM. Each of these sections can be characterized by minimal diameter of the 
ellipsis (c?min)fc and the angle (<po) k of its rotation in the cross-sectional plane, while the 
maximum diameter rf max of the ellipse is the same for the whole sample. The values of 
these parameters are found from experimental set of diameters {{di) k } minimizing the mean 
squared deviations 

with respect to d max and the set of N individual diameters (d min ) k and angles (y?o)fc for 
each of N sections. Here d\ (cp) is given by Eq. (J7J) with corresponding parameters (d min ) k 
and (<po) k - The solution of these minimum conditions is quite cumbersome and rendered in 
Appendix. Here we show the result 



d 2 

"max 



4 . 2 f— — 1 



« in ) k = -A k - - y/Bl + CI - -d 2 max , (20) 

(<Po) k = ^arctan^, (21) 
where A k ,B k and C k are Fourier transforms of squared diameters of the specimens: 

•4- 5 (22) 

^'iE",^™^.). (23) 

2 M 

^^Lt^w- (24) 



IV. INTERPRETATION OF DATA 



In this section we present observable dependences of geometrical sizes of samples on 
fluence $. In the following we show only curves averaged over the row of specimens with the 
same orientation (|| and _L), diameter (6 and 8 mm) and average irradiation temperature 
T. Samples of assembly A were exposed to radiation in reactor Bor-60 at the temperature 
T = 460 ± 25°C, samples of assemblies B and C have irradiation temperatures 550 ± 25°C 
and 640 ± 25°C, respectively. 
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a) Length changes. Most information about radiation-induced shape-change of 
graphite is usually obtained using data describing the change of length L of samples with 
neutron fluence $. For example, the volume change AV is commonly estimated by the for- 
mula ((H), using the data on length variations for specimens of parallel, (AL) , and perpen- 
dicular, (AL)" 1 , orientations. The dependence of relative elongations (AL/L)" and (AL/L) ± 
of studied samples on fluence $ is shown in Fig. [3j As one can see, the change of the sample 




FIG. 3: The dependence of relative length changes (in %) for samples (ad) from the assembly 
a = A,B,C with diameter d (6 mm - thin lines, 8 mm - thick lines) with parallel (||) a) and 
perpendicular (_L) b) orientation oj on neutron fluence <3? (E n > 0.18 Mev). 

size strongly depends on sample shape because of different constraints imposed on domains 
changing their shape in samples of different diameter d and orientation u. The radiation 
induced shape-change of domains is only weakly restricted in specimens with perpendicular 
orientation and of smallest diameter d = 6 mm, while the domains experience strongest 
restrictions in specimens with parallel orientation and of diameter d = 8 mm. 

b ) Volume changes. In Fig. H]we show the dependence of the sample volume on fluence 

$. 

The evolution of the sample volume with the fluence $ can be calculated in two different 
ways: using variations of macroscopic sample sizes (Eq. (fT6l) ) and using the change of local 
sizes of domains (Eq. f|T7|) ). As one can see from figure 0] a), both ways give essentially the 



11 




FIG. 4: Relative volume change AV/V (in %) as function of fluence $ for assembly A with 
temperature T = 460 ± 25°C. a) solid lines (-F) show (AZ/Z)|| + 2 (A//Z)j_, dash-dotted lines show 
AL/L + 2 (Ad/d) for 6 mm samples of both parallel (||) and perpendicular (_L) orientations. Since 
similar curves are obtained for all other samples, in figures b) and c) we show only single dependence 
of AV/V = AL/L+2 (Ad/d) on <3? for assemblies A, B and C (solid lines, with average temperatures 
460,550 and 640°C, respectively), and also asymptotic dependence of these curves C[(Al/l) ± — 
(Al/l)^] 2 (dashed lines, C ~ 7 for all curves) at large fluence <3?. 

same result: 

(AV/V) F = AL/L + 2 (Ad/d) = (A///),, + 2 (Al/l) ± (25) 

We observe that the volume of the specimen only slightly varies with its orientation u. 
Therefore, values of principal strains (AZ/Z)|| and (Al/l) ± , through which the volume change 
is expressed, can be considered as universal characteristic of graphite, weakly depending on 
the specimen orientation u. In contrast to the volume change (see Fig. H]) the length L 
and the diameter d of the sample (see Fig. |3]) are not universal and strongly depend on the 
sample shape. These values depend not only on principal strains, but also on the flattening 
factor e which is not universal because of different constraints imposed on domains changing 
their shape under irradiation in small specimens of different orientations. We observe that 
variations of the volume are depressed for samples of larger diameter, d = 8 mm, because of 
stronger constraints imposed on shape-changing domains in larger samples. 

The dependence of volume change on fluence $ is strongly non-monotonic: at small 
fluences $ < $ m i n the graphite shrinks with dose $, and it dilates at $ > $ m i n . Dotted lines 
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in Figs. H] a) and b) show the dependence 



C (A///) 



(A///) 



(26) 



with coefficient C = 7. As one can see, this dependence describes relatively well asymptotic 
behavior at large $ > $ m i n of volume changes in graphite for all temperatures, diameters 
and orientations. The knowledge of the asymptotic behavior is important both to under- 
stand undergoing mechanisms of graphite deformation under irradiation and to extrapolate 
experimental data in a region of higher radiation dozes not reachable for current experiments. 

Using Eq. ( Fl2|) we find the following expression for combined equation (JTJ) traditionally 
used to describe the volume change AV C during irradiation: 



We observe, that Eq. (127|) gives different results (AV C < AV F ) from real relative change 
of sample volume given by Eq. ( |25|) (see dotted line in figure BJ due to ellipticity, e > 0, of 
specimen cross-sections. We conclude that the traditional expression ( 1271) can not be used 
to calculate the actual change of volume of graphite samples. The replacement of actual 
volume change AV F by AV C can be safety applied only in the case of homogeneous texture 
of graphite, see Fig. [1] a. Disorder in domain orientations (Fig. [1] b) leads to significant 
difference between commonly used (Eq. f )27j) ) and actual (Eq. (125]) ) volumes of samples. 

c) Parameters of domains. The orientation of domain anisotropy axis is determined 
by two angles: the azimuth angle (po and the polar angle 6, see Fig. [2]b). Below we estimate 
parameters of domains from the knowledge of statistics of these angles for samples under 
consideration. 

Azimuth angle <£>o determines orientation of the axial axis in the cross-sectional plane. 
The angles {<fo) k are random variables of the cross-section number k, that are correlated at 
the distance about the domain size £. In general, the amplitude of angle fluctuations can 
be different for different samples. In order to exclude the effect of such amplitude variations 
consider the ratio 

2r 

r =k (28) 

of two sums for each of samples 



(AV/Vf = (AV/V) F - {2e L - e 11 ) < (AV/V) F 



(27) 




N 



N 




(29) 



k=l j=k+l 



k=l 
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The first sum is going over all pairs of cross-sections (an analog of a gyration radius), while 
the second one is only over nearest pairs (an analog of end-to-end distance of a polymer 
chain). Since each of these terms is quadratic form of angles (<po) k , their ratio (f28l) does not 
depend on amplitude of fluctuations and encrypts only information about domain size. If 
there are no correlations between angles fk this ratio tends to 1 in the limit of large N — > oo, 
when we can change [((p ) k - (v^o)^ 2 in Eqs. ([29]) by its average, [{<po) k - (y^] 2 = In 
order to determine the dependence of the parameter r on domain size £, consider the case 
when the group of I consecutive angles f n i+j, 3 — 1, • • • , / has the same value <p n for integer 
n, I, whereas the values tp n for each group n are independent random variables. At large 
N — > oo changing [(ipo) k — (v 9 o) J ] 2 in Eq. (I2"9~j) by its average (0 for k and j belonging to 
the same group and 2ip^ otherwise) we find that the parameter r equals the length I of the 
group. Therefore, we can identify £ = rL Qk with the average domain size along the direction 
of the cylinder axis, L ok is the distance between different sections of the sample, see Fig. |2j 
In Table U we present the results of calculations of the parameter r for specimens of parallel 



TABLE I: 



Assembly (°C) 


A (460) 


B (550) 


C (640) 


Diameter (mm) 


6 


8 


6 


8 


6 


8 


r" 


1.2 


1.2 


1.3 


1.1 


1.1 


1.4 


r ± 


1.6 


1.4 


1.4 


1.3 


1.3 


1.5 



(r") and perpendicular (r -1 ) orientations averaged over all such samples. 

Since the distance between different cross-sections of the specimen is L fe = 4.5 mm, we 
get estimation £ ~ 6 mm for the typical domain size. As shown in Table [H specimens of 
perpendicular orientation have higher r than specimens of parallel orientation (r 1 - > r"). The 
larger longitudinal size of domains with perpendicular orientation is related to less restricted 
environment of such domains, that weaker constraints the radiation-induced evolution of 
domain shape. 

The orientation of local axial axis with respect to the global axis of the cylinder is 
described by the polar angle 6 (see Fig. |2] c). Combining Eqs. ( fl3|) and ( [181) we find 
the average 
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Analyzing experimental data we show in Fig. [5] the dependence of sin 2 9 on fTuence $ for 
samples of different diameters d — 6,8 mm and orientations u (|| and 1). 



sin u y — on \ 

o- 1.0- 




k - : — 277 / fill 

sin e c / 6 
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81 
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0.6- 
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^81 . 2 . 


61 

81 


81 

1 - ^> 



1. OxlO 22 2 - 1. X1O 22 2. o l. (J, xl0 22 2. 



FIG. 5: The dependence of sin 2 9 on fluence $ for samples of assemblies a = A, B,C with different 
average temperatures. Symbols duj correspond to samples of diameter d = 6 and 8 mm and 
orientation uj (|| and 1). 

One may ask whether the axial axis can be rotated as the result of radiation induced 
shape-change of domains? Inspecting Fig. [5] we conclude that this is not the case, and the 
orientation angle 9 does not vary with the fluence $. The statistics of this angle only weakly 
depends on temperature T during irradiation, and thus is mainly determined by fabrication 
conditions. An important special case are samples of diameter d = 6 mm with perpendicu- 
lar orientation, that are strongly inhomogeneous and have regions where radiation-induced 
shape-change is only weakly restricted by the environment. Such heterogeneity of the sam- 
ple is the reason of apparent rise of sin 2 9 above 1 with the irradiation dose in Fig. [5j This 
effect is absent for samples of diameter d = 8 mm, supporting our previous estimation £ ~ 6 
mm for the domain size. 

d) Principal strains. The dependencies of principal strains (AZ/Z)|| and (Al/l) ± on 
neutron fluence $ for all studied samples are collected in Fig. |6j Below we use these data 
to discuss the dependence of these strains on shape and size of samples, while the dose 
dependence will be considered later. 

The dimension effect is the most pronounced for parallel principal strain (Al/l),,, see 
Fig. |6]a). At small fluence $ the domain always shrinks in parallel direction. The amplitude 
of the shrinking (Ai/Z),, is maximal for specimens with perpendicular orientation, since 
domains of such samples are less restricted to grow with increasing dose $ with respect to 
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FIG. 6: The dependence of principal strains (AZ/Z)ii and (A///)j_ (in %) on neutron fluence for 
samples (adco) in the assembly a = A,B,C with diameter d (6 mm - thin lines, 8 mm - thick lines) 
and orientation oj (|| - solid lines, _L - dash lines). The same notations as in Fig. 

specimens of parallel orientation. This effect strongly depends on the sample size. While 
the strain (A///)|| dramatically changes with orientation u for samples of diameter d = 6 
mm, the value (Al/l),, relatively weakly changes for samples of diameter d = 8 mm. This 
conclusion is in agreement with obtained above estimation of the domain size £ ~ 6 mm, 
and we expect that the dimension effect is small only for samples of diameter d > 8 mm. 

The variation of (Al/l)^ with the fluence $ is related to corresponding growth of microc- 
racks which is hampered by internal stresses in graphite. Such constrained effect of internal 
stresses is minimal for specimens with perpendicular orientation. Similar behavior is ob- 
served for perpendicular principal strain (Al/l) ± , as well. At small $ the value of (Al/l) ± 
weakly depends on fluence The domain shrinks in perpendicular direction at small tem- 
perature T < 600° C and extends at high temperatures T > 600° C. The principal strain 
(Al/l) ± strongly varies with orientation u only for samples of diameter d = 6 mm, and do 
not depend on orientation for samples of diameter d = 8 mm. The amplitude of variation 
of (Al/l) ± with neutron fluence $ is maximal for the least constrained in perpendicular 
principal direction specimens with parallel orientation, and is minimal for specimens with 
perpendicular orientation. 

The amplitude of anisotropy increases with the rise of the irradiation temperature T. 
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This effect is related to the development of an ensemble of anisotropic microcracks and the 
rise of the microcrack volume with increasing temperature T. The increase of the volume 
of microcracks with the rise of the temperature is responsible for corresponding rise of 
the sample volume, clearly observed in Fig. HI At high temperature (T = 650° C) the 
contribution of microcracks is so big, that at low doses it prevents the shrinkage of the 
specimen length in parallel principal direction, (Ai/i),,, although there is always a region 
of initial volume shrinking since at small $ the domain shrinks in perpendicular principal 
direction, (Al/l) ± . 

Principal strains (Al/l),, and (Al/l) ± demonstrate essentially non-monotonic dependence 
on fluence $, see Fig. [6j In order to understand the reason of such behavior we can present 
these strains in the form 



Al 
T 



Al\ 



1AV 



3 V 3 
1 AV 1 

TV + 3 



Al 

T 

Al 

T 



AZ 

T 

Al 

T 



(31) 
(32) 



The non-monotonic dependence of these strains on $ is related to corresponding non- 
monotonic dependence of the volume, see Fig. 0J The dependence of the difference of 
principal strains, (Al/l) ± — (Al/l)^, on fluence $ is monotonic and can be extrapolated 
by quadratic function, see Fig. a: 



(AZ/Z) ± - (AZ/Z),, ~ 



(33) 



The shape of samples becomes more and more crankshaft- like with the rise of the neutron 
fluence The amplitude of this effect is characterized by the average flattening factor e of 
the elliptic cross-sections of the specimen, Eq. ( fl"3l) . Inspecting experimental data we find 
that e grows quadratically with the fluence $: 



£($) ~&$ 2 



(34) 



The dependence (13"1|) is illustrated in Fig. b), where we plot e for different assembles, 
diameters and orientations of specimens as function of $ 2 . 

The coefficients Pi and /3 e for different orientations u, diameters d and assembles a are 
shown in Table IIVI The factor $ monotonically grows with the temperature T during 
irradiation (from assembly A to C). The difference of fix for specimens of parallel, ($)", 
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FIG. 7: The dependence (in %) of the principal strain difference (Al/l) ± — (AZ/Z).. a) and of the 
flattening factor e b) on squared fluence $ for all samples. The same notations as in Fig. [3l 
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and perpendicular, (/?;) , orientations decreases with the rise of the sample diameter, but 
the gap between these values exists even for maximal diameter d = 8 mm. Therefore, we 
conclude that in our experiment even the most thick samples have elastic properties different 
from properties of a bulk graphite. 

The coefficients (3 E increase with the rise of the temperature T. According to Eq. ( TT3]) 
the two coefficients 0i and (3 £ are proportional to each other: 

& = tifeft (35) 

Higher values of (3 e for specimens of perpendicular orientation with u = 7r/2 are related 
to larger values sin 2 9 with respect to the case of parallel orientation and orientation angle 
u = 0. We expect that for bulk graphite both coefficients vanish, (3i = (3 E = 0, and there 
should not be any difference between the two ways of measurement of the volume change, 
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AV C = AV F , Eqs. (ESJ and (J2TJ) . 

At first sight simple quadratic dependence f )33|) looks very surprising, taking into account 
complex non-monotonic dependence of volume changes on the fluence $, see Fig. HJ at low 
dose $ < $ m in the sample is shrinking, while at high doses it dilates. It is usually assumed 
that appearance at large fluences $ > $ m i n of growing branches in volume and elongation de- 
pendences is related to formation of new set of microcracks in graphite at moderate fluences 
$ > $mi n .— Simple dependence f )33|) demonstrates that really nothing special is happened 
with graphite above the cross-over $ m j n . In Ref.— we proposed alternative explanation of 
such non-monotonic dependences as the result of the interplay of radiation induced changes 
of the shape and variations of elastic moduli of crystallites in polycrystal graphite. At low 
doses, $ < $ m i n , the change of the shape of domains wins, leading to initial shrinking of 
graphite samples. At high doses, $ > $ m i n , the main contribution to graphite deformation 
comes from the variation of elastic moduli, leading to dilation of graphite. 

e) Relative contributions. In Fig. [8] we show the relative contribution of four impor- 
tant geometrical characteristics of samples of assembly A: the change of length L, the change 
of diameter d and two principal strains (Al/l)» and (Al/l) ± as functions of fluence Due 
to Eq. fl25|) there is linear dependence between these variables. We observe that measured 
elongations L and d lie in between principal elongations (A///)|| and (Al/l) ± because of 
constraints imposed on the domain shape in the volume of graphite. 

V. THEORY: INTERNAL STRESSES IN GRAPHITE 

Intuitively, when all domains are strictly oriented, the radiation-induced change of their 
shape does not cause any internal stresses. According to the terminology of Ref ~ such 
domains are in the state of a free dilatation. In this state the relative volume change 
of a sample of arbitrary shape and of the macro-graphite are the same. This idealized 
case corresponds to expression ([I]) and the upper curves in Figs. in the sense that they 
describe the change of domain shape in almost no stress, despite of apparent disorientation of 
domains. In this case internal stresses are absent because the size of domains is comparable 
with the minimum size (the diameter) of the samples. In the real case of macro-graphite with 
all dimensions much larger than the domain size, the situation is fundamentally different. 
Below, we consider this case and describe the emerging internal stresses because of domain 
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A 6 mm, _L 



15- 




A 8 mm, _L 



FIG. 8: The dependence of relative elongations (in %) of samples in reactor BOR-60: AL/L 
(dashed thin lines, symbol L), Ad/d (solid thin lines, symbol d), (A///)j_ (solid thick lines, _L) and 
(Al/l),, (dashed thick lines, |[) on fluence <J? for samples of the assembly A with average irradiation 
temperature T = 450° C. The same notations as in Fig. [31 



disorientation (see Fig. [T] b) following the method developed in Ref .— . 

Let us analyze the meaning of differences in the volume change (AV/V) obtained in 
Eqs. ( Il7p and ([T]). Expression ffTTj) corresponds, as already mentioned, to the radiation- 
induced "free dilatation". Now turn to Eq. (JT]). Longitudinal elongation in this expression 
(for example, (AL/L),, for parallel orientation) describes the real strain of the sample, 
taking into account the disorientation of its domains. Transversal elongation ((AL/L) ± 
for parallel orientation) in expression (DQ) corresponds to totally suppressed deformations 



20 




FIG. 9: The dependence of relative volume changes (in %) for samples of diameters d = 6 and 8 mm 
on fluence $ for assemblies a = A,B, C, corresponding to irradiation temperatures T = 460, 550 
and 640°C: upper dash-dotted lines F give actual values, (AV/V) , Eq. (fl~6j) , lower dotted lines 
C correspond to traditional expressions, (AV/V) C , Eq. CO), and middle solid lines M are our 
predictions for bulk graphite, (AV/V) , Eq. ([38]) . 

of adjacent domains of distinct orientations. Eq. ([1]) gives the relative volume change of a 
simply connected macro-graphite assembled from the samples of say, parallel orientation, 
when all discrepancies between adjacent shape- changing domains are straighten by stresses 
at domain boundaries. The difference in values (ITTj) and determines the excess volume 
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fraction of such overlapping radiation-induced regions: 

c fAV\ fAV\ F fAV\ c n , II 

H^H^J "(IT) (36) 
After reaching the equilibrium at internal boundaries between the domains the excess volume 
fraction 5q(AV/V) partially relaxes. Considering a macro-graphite as a body, containing an 
ensemble of equilibrium dilatating objects whose dimensions are much smaller than the size 
of the body with free boundaries, one can find the dilatation component of the total strain. 
Following Ref.— we get in the leading order in (AV/V): 

Tr{u) ^1^(2^-5-") (37) 

Since in the reference state the relative volume change is (AV/V) C , we find expression for 
the relative change of the volume of macro-graphite: 

{AV/V) M = (AV/V f + Tr («) 

/ AL\ " / AL\ ^ 21 -2v , | || \ , , 

+2 {-) + 3— ( & - £HI ) (38) 

The relation between different definitions of relative volume change, (AV/V) F , (AV/V) C 
and (AV/V) M , is shown in Figs. [9] taking into account the value of Poisson's ratio v = 0.2 
for graphite GR-280. 

The obtained result is in agreement with frequently noted in the literature experimental 
fact - that large samples of extruded graphite change their shape stronger than the thinner 
samples. For example, the paper— presents the results of measurements of irradiated sam- 
ples, ~ 10 cm x 10 cm x 61 cm, and also of cylindric specimens with diameter 1.09 cm 
and length 10 cm, in the case of parallel and perpendicular orientations with respect to the 
extrusion axis. The observed values of the volume change (the shrinkage) for larger sam- 
ples are approximately twice as large. Using expressions fl36|) and (138]) the relation between 
macro-graphite and thin samples can be written as 

VI. CONCLUSIONS 

In this work we present data on both longitudinal and transverse measurements of 
radiation-induced shape-change for graphite specimens of cylindrical shape. These data 
show that: 
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• The cross-sections of the specimens change their shape under irradiation from initial 
circular to elliptical one, and orientations of main ellipse axes randomly vary along 
the sample. 

• The relative volume changes, obtained by direct measurement of AV/V F and using the 
traditional expression (AV/V) in terms of relative length variations of samples with 
parallel and perpendicular orientations, are essentially different (see expression (1271) 
and Fig. E]). 

To explain these facts, authors proposed a new model of the morphology of extruded 
graphite, considering it as an ensemble of dilated domains with orientations disoriented 
with respect to the global axis of symmetry (direction of extrusion) of macro- graphite. In 
this model the symmetry properties of both macro-graphite and domains are the same - the 
transversal isotropy. 

For each cross-section of cylindrical specimens of diameters 6 and 8 mm with parallel and 
perpendicular orientations we present results of measured values (AL/L)" , (AL/L) ± , <i max 
and d min . We derive closed system of equations for parameters (AZ/Z),, , (Al/l) ± and sin 2 9 
characterizing deformation of domains and their orientation and calculate these values. 

Our model predicts that domains of samples with transverse dimensions on the order of 
domain size should experience free (unconstrained) radiation-induced dilation. In this case, 
the values (AZ/Z)y , (Al/l) ± and sin 2 # must match for samples with parallel and perpendic- 
ular orientations. This condition is satisfied for specimens of diameter 6 mm and 8 mm (see 
Fig. [8]). The data for such free-dilatating domains were used in the subsequent analysis and 
evaluations. 

The main results can be summarized as: 

a. In macro-graphite radiation-induced dilation of disoriented (sin 2 9 > 0, see Fig. [T]b) 
domains leads to the development of cross-domain stresses. Of particular interest is 
the dilatation component of the stress leading to the difference between real volume 

F C 

change (AV/V) and commonly used expression (AV/V) for this change. Note that 
in the case of completely ordered domains (see Fig. [1] a) internal stresses disappear 
and all expressions for the volume change are the same, 




(40) 
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(3. It is shown that the shrinkage curves obtained both by direct measurements and by 
calculation of relative size changes do not describe the real shrinkage curve of macro- 
graphite. We derived the relation 

— ) =(-) +2 {-) + s—^-^ (41 > 

that allows to obtain the shrinkage curve of macro-graphite by direct measurements 
of strains in samples of small sizes. 

7. Our analysis leads to important conclusion about the presence of size restrictions on 
experimental specimens. In the case when minimum sample -^miii IS 

much larger 

than the domain size £, L min ^> £, the elongation measurements directly match the 
behavior of macroscopic graphite. However, such a condition is very hard to implement 
in reactor experiments. 

The simplest realized condition is to satisfy the ratio L min ~ £, with subsequent conversion 
of obtained data using the relation (1411) . In our experiments, as seen from Fig. [SJ, the 
condition L min ~ £ is realized for the set of specimens with diameters of 6 mm and 8 
mm. This criterion is in a satisfactory agreement between principal elongations (Al/l),< and 
(Al/l) ± of dilatating domains of different orientations. 

Acknowledgments 

The authors consider as their pleasant duty to thank B.N. Bennesch, A.N. Maltsev and 
A.V. Nechitailo who made a significant contribution to the experiments. 



Electronic address: subbotin@atomtechnoprom.ru 



' Electronic address: panyukov@lpi.ru 



* Electronic address: ivanov@lpi.ru 



1 John H.W. Simmons, Radiation damage in graphite, Oxford, Pergamon Press (1965). 

2 B. Kelly, Radiation damage in graphite and its relevance to reactor design, Progress in Nucl. 
Energy, 2 (1978) 219-269. 



24 



3 Irradiation damage in graphite due to fast neutrons in fission and fusion systems, IAEA- 
TECDOC-1154, September (2000). 

4 G. Haag, Properties of ATR-2E Graphite and property changes due to fast neutron irradiation, 
Forschungszentrum, Jiilich, (2005). 

5 Ya.I. Shtrombakh, B.A. Gurovich, P.A. Platonov and V.M. Alekseev, Radiation damage of 
graphite and carbon-graphite materials, J. of Nucl. Mater. 225 (1995) 273-301. 

6 A.R. Ubbelohde and F.A. Lewis, Graphite and its crystal compounds, Oxford, Clarendon Press, 
1960. 

7 W. Windes, T. Burchell, R. Bratton, Graphite technology development plan, INL/EXT-07- 
13165 (2007). 

8 J.F. Nye, Physical properties of crystals, Oxford, Claredon Press, (1964). 

9 A.E.H. Love, Mathematical theory of elasticity, London, (1952). 

A.V. Subbotin, O.V. Ivanov, I.M. Dremin, V.P. Shevel'ko, V.A. Nechitailo, Primary radiation 
damage to graphite in different types of reactors. Atomic Energy, 100 (2006) 199-209. 

1 S.V. Panyukov and A.V. Subbotin, Atomic Energy, 105 (2008) 32-35. 

2 J.D. Eshelby, Solid State Physics 3, pp. 79-144. Academic Press, New York (1956). 

3 R. Nightingale and E. Woodruff, Radiation-Induced dimensional changes in large graphite bars, 
Nucl. Sci. Eng., 19 (1964) 390. 

Appendix 



The angular dependence of the diameter ([7]) of deformed cylinder with initial circular 
cross-section can be rewritten in the form 



Solution of minimum conditions 



d 2 (<f) = a k + b k cos (2ip) + c k sin (2ip) 



(Al) 



where 




(A2) 



2 [^max - Klin) J ^ [ 2 {<fo) k } 
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Knowing coefficients a k ,b k and c k we can calculate from these equations the two principal 
diameters d max and (d min ) k : 

1 N 

rf Lx = T7 ^ H + ^ cos ' 2 + c k sin [2 (tpo) k }} , (A3) 



N 

k=l 



Kin) fc = 2a fc -d: 



2 

max 



In order to find coefficients a k , b k and c k from experimental data we substitute Eq. (lAlj) 
into Eq. ([7j). The result of this substitution can be simplified using normalization conditions 

M M 

sin 2 (2caA = ■ 



cos 2 (2^) = sin2 ( 2 ^) = T (A4) 

i=l i=l 



and orthogonality relations 

M M M 

cos (2<£?j) sin (2y?i) = cos (2^) = sin (2</?j) = (A5) 
After some algebra we get 

N 

G 2 = MY^ [2a 2 +b 2 k + c 2 k - 4a k A k - 2b k B k - 2c k C k ] , (A6) 

k=l 

where A k , B k and C k are defined in Eq. (|T9|) - (12TI) . Minimizing expression (1A6j) with respect 
to (d min ) k , d max and ((po) k we find equations for these variables 

2 (a fe - A k ) - (b k - B k ) cos [2 {<f ) k ) - (c k - Ck) sin [2 (<fo) k ) = 0, 

EN 
t {(a k - A k ) + (b k - B k ) cos [2 {<p ) k ) + (c k -C k ) sm[2( ( p ) k }} = (A7) 

- (fo fc - B k ) sin [2 (</? ) fc ] + (cfe - C fe ) cos [2 {<fo) k ) = 
Inspection of the above equations shows that b k and c k can be written in the form 



h = \jb 2 k + c\ cos [2 (p„) J , c fc = y'o 2 + c 2 sin [2 (p„) J (A8) 
where (<po) k is the solution of equation 

B k sin [2 (^q) J = C k cos [2 (^ ) J (A9) 
Substituting Eqs. (jSSJ) and (TA9]) into Eq. (J2Ij) we find 



2 (a* - A fe ) - y'o 2 + c 2 + ^ 2 + C\ = 0, (A10) 

2 EL - ^ + - v^+^} = (A11) 



Combining these equations we get 
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hi + c = Bl + C 



(A12) 



and also 



a k = A k + X /Bl + CI 



Substituting Eq. ( 1A12I) into Eq. ( 1A8I) we find coefficients b k and c^: 



(A13) 



ki ' 



^ + C fe 2 sin[2 fa,), 



(A14) 



As the last step of our derivation we substitute Eqs. (IA13I) and (1A14I) into Eqs. (1A3[) and 
find diameters <i max and (d min ) k , Eqs. (TT^|) and (I2"0j) . 



